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Lattice QCD with domain wall quarks and applications to weak matrix 
elements 

T. Blum* and A. Sonit a 

a Department of Physics, Brookhaven National Laboratory, Upton, NY 11973 

Using domain wall fermions, we estimate Bk{/J. ~ 2 GeV) = 0.628(47) in quenched QCD which is consistent 
with previous calculations. At 6/g 2 = 6.0 and 5.85 we find the ratio fx/wip in agreement with the experimental 
value, within errors. These results support expectations that 0(a) errors are exponentially suppressed in low 
energy (E <C a~ ) observables, and indicate that domain wall fermions have good scaling behavior at relatively 
strong couplings. We also demonstrate that the axial current numerically satisfies the lattice analog of the usual 
continuum axial Ward identity. 
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A basic feature of the strong interactions 
has been missing in lattice calculations, the 
SU(Nf )l x SU(Nf )r chiral flavor symmetry of the 
light quarks. We recently reported Q on calcula- 
tions using a new discretization for simulations of 
QCD, domain wall fermions (DWF) §,§, which 
preserve chiral symmetry on the lattice in the 
limit of an infinite extra 5th dimension. There it 
was demonstrated that DWF exhibit remarkable 
chiral behavior Q even at relatively large lattice 
spacing and modest extent of the fifth dimension. 
Here we give further results using DWF which are 
of direct phenomenological interest 

In addition to retaining chiral symmetry, DWF 
are also "improved" in another important way. 
In the limit that the number of sites in the extra 
dimension, N s , goes to infinity, the leading dis- 
cretization error in the effective four dimensional 
action for the light degrees of freedom goes like 
0(a 2 ). This theoretical dependence is deduced 
from the fact that the only operators available to 
cancel 0(a) errors in the effective action are not 
chirally symmetric. For finite N s , 0(a) correc- 
tions are expected to be exponentially suppressed 
with the size of the extra fifth dimension. Our 
calculations for Bk show a weak dependence on 
a that is easily fit to an a 2 ansatz. 

We use the boundary fermion variant of DWF 
developed by Shamir. For details, consult Ka- 



plan Q and Shamir See Ref. || for a discus- 
sion of the Ad chiral Ward identities (CWI) sat- 
isfied by DWF. Our simulation parameters are 
summarized in Table |l|. 

Table 1 

Summary of simulation parameters. M is the five 
dimensional Dirac fermion mass, and m is the 
coupling between layers s = and N s — 1. 





size 


M 


m 


# conf 


5.85 


16 a x 32 x 14 


1.7 


0.075 
0.05 


34 
24 


6.0 


16 3 x 32 x 10 


1.7 


0.075 
0.05 
0.025 


36 
39 
17 


6.3 


24 a x 60 x 10 


1.5 


0.075 
0.05 


11 

14 



We begin with the numerical investigation of 
the lattice PCAC relation. The CWI are satis- 
fied exactly on any configuration since they are 
derived from the corresponding operator identity. 
We checked this explicitly in our simulations. In 
the asymptotic large time limit, we find for the 
usual PCAC relation 
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which goes over to the continuum relation for 
am, <C 1 and iV s — > oo (see Ref. || for opera- 
tor definitions). The second term on the r.h.s. is 
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anomalous and vanishes as N s — > oo. It is a mea- 
sure of explicit chiral symmetry breaking induced 
by the finite 5th dimension. At 6/g 2 = 6.0 and 
N s = 10 we find the l.h.s. of Eq.|to be 0.1578(2) 
and 0.1083(3) for m = 0.075 and 0.05, respec- 
tively. The anomalous contributions for these 
two masses are 2x(0.00385(5) and 0.00408(12)), 
which appears to be roughly constant with m. In- 
creasing N s to 14 at m = 0.05, the anomalous 
contribution falls to (2x) 0.00152(8) while the 
l.h.s. is 0.1026(6), which shows that increasing 
N s really does take us towards the chiral limit. 

Next we investigate the matrix element of Oll- 
At 6/g 2 — 5.85 the two data points extrapolate 
linearly to -0.0005(100) at m = 0. At 6/g 2 = 6.0 
the three data points extrapolate to -.004(9). At 
6/g 2 = 6.3, the two points extrapolate to 0.05(3). 
This slight overshoot is not unexpected since the 
values of m used here correspond to rather heavy 
quarks. In our initial study we found a similar be- 
havior 0], and as the quark mass was lowered, the 
required linear behavior set in. All of the above 
results are for N s = 10 except at 6/g 2 = 5.85 
where N s = 14 was used for reasons explained 
below. 
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Figure 1. The kaon B parameter. The cross (not 
used in the fit) denotes the partially unquenched 
result using KS sea quarks ||. 

In Fig. H] we show the kaon B parameter. The 



results for Bk depend weakly on 6/g 2 , and are 
well fit to a pure quadratic in a. We find Bk([J. = 
a" 1 ) = 0.628(47) in the continuum limit. This 
value is already consistent with previous results j7|] 
though it does not include the perturbative run- 
ning of Bk to a common energy scale. This re- 
quires a perturbative calculation to determine the 
renormalization of Oll, which has not yet been 
done. The energy scale at 6/g 2 = 6.0 is roughly 
2 GeV. Also, we find B K (ji = a" 1 ) = 0.67(4) on 
a set of 20 Kogut-Susskind lattices with ttiks = 
0.01 and 6/g 2 — 5.7|| and the same five dimen- 
sional lattice volume as the point at 6/g 2 = 6.0. 

At 6/g 2 = 6.0, we have also calculated Bk us- 
ing the partially conserved axial current A"(x) 
(and the analogous vector current). This point 
split conserved current requires explicit factors of 
the gauge links to be gauge invariant. Alterna- 
tively a gauge non-invariant operator may be de- 
fined by omitting the links; the two definitions 
become equivalent in the continuum limit. Re- 
sults for the gauge non-invariant operators agree 
within small statistical errors with those obtained 
with naive currents, Fig. [j](see Ref.||,|l[ for oper- 
ator definitions). The results for the gauge in- 
variant operators are somewhat larger: B K lv (/.i = 
a" 1 ) = 0.857(20) and 0.946(28) at m = 0.05 and 
0.075, respectively. A similar situation holds in 
the Kogut-Susskind case where it was shown that 
the gauge invariant operators receive appreciable 
perturbative corrections which bring the two re- 
sults into agreement [||. 

Assuming PCAC, the pseudoscalar decay con- 
stant is determined from the measurement of 
(0|J 5 a |P). We find f K = 159(14) MeV and 
164(12) MeV for 6/g 2 = 6.0 and 5.85, respec- 
tively. The errors are statistical and do not in- 
clude the error in the lattice spacing determina- 
tion from am p . The central values agree with 
experiment, = 160 MeV. The lattice spac- 
ing determinations from am p are a -1 = 1.53(27), 
2.09(21), and 3.20(81) GeV at 6/g 2 = 5.85, 6.0, 
and 6.3, respectively. Alternatively, we may form 
the dimensionless ratio /k/to p . We find for 
6/g 2 = 5.85 and 6.0, f K /m p = 0.213(42) and 
0.206(27), where we have added the statistical 
errors naively in quadrature. The experimental 
result is 0.208. The decay constant calculated di- 
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rectly from the matrix element of the partially 
conserved axial current agrees with results us- 
ing the matrix element of the pseudoscalar den- 
sity, up to the anomalous contribution which is 
slightly less than the statistical error. While the 
above indicate good scaling behavior, they must 
be checked further with improved statistics and 
a fully covariant fitting procedure. More impor- 
tantly, the continuum limit still has to be taken: 
a recent precise calculation using quenched Wil- 
son quarks by the CP-PACS collaboration gives 
a value for /if/m p in the continuum limit that 
is inconsistent with experiment ]l0[ ] , so the above 
agreement with experiment may be fortuitous. 
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Figure 2. The pion mass squared. 

In Fig. |^ we show the pion mass squared as 
a function of m. The data at Q/g 2 — 6.0 and 
6.3 are consistent with chiral perturbation the- 
ory. At 6/g 2 = 5.85 the two masses extrapo- 
late to 0.045(10) for N s = 10 and 0.031(13) for 
N s = 14. This discrepancy is probably not due to 
higher order terms in the chiral expansion since 
the physical quark masses are light compared to 
the masses at the other couplings, and the curva- 
ture would have the wrong sign. We see a large 
downward shift in m 2 as N s goes from 10 to 14. 
However, increasing N s to 18 at m = 0.075 has 
a negligible effect(the point at m = 0.05 should 



also be checked since the slope may increase) . We 
note that the anomalous contribution to Eq. |l| is 
more than double the value at 6/g 2 = 6.0. In the 
case of the vector Schwinger model, it was found 
that topology changing gauge configurations can 
induce significant explicit chiral symmetry break- 
ing effects Further investigation is required. 

Our study shows that DWF are an attractive 
alternative for lattice QCD calculations where 
chiral symmetry is crucial. DWF yield good 
agreement with expectations from chiral pertur- 
bation theory without the complicated mixing of 
operators required with Wilson quarks, or the 
entanglement of flavor and space-time degrees 
of freedom as with Kogut-Susskind quarks. Up 
to exponentially small corrections, DWF main- 
tain the full chiral symmetry of QCD at rel- 
atively strong couplings, and thus should have 
more continuum-like behavior. The data pre- 
sented here seem to indicate just that, though 
future studies with improved statistics are needed 
to confirm this. This improved scaling may com- 
pensate for the added cost of the extra dimension. 
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